Exploring the dark accelerator HESS J1745-303 with Fermi Large 

Area Telescope 



C. Y. Hui\ E. M. H. Wu^, J. H. K. Wu^, R. H. H. Huang^, K. S. Cheng^, P. H. T. Tam^, 

and A. K. H. Kong^'^ 



ABSTRACT 

We present a detailed analysis of the 7— ray emission from HESS J1745-303 with the data 
obtained by the Fermi Gamma-ray Space Telescope in the first ^ 29 months observation. The 
source can be clearly detected at the level of ~ 18a and ^ 6cr in 1 — 20 GeV and 10 — 20 GeV 
respectively. Different from the results obtained by the Compton Gamma-ray Observatory, we 
do not find any evidence of variability. Most of emission in 10 — 20 GeV is found to coincide with 
the region C of HESS J1745-303. A simple power-law is sufficient to describe the GeV spectrum 
with a photon index of F 2.6. The power- law spectrum inferred in the GeV regime can be 
connected to that of a particular spatial component of HESS J1745-303 in 1 — 10 TeV without any 
spectral break. These properties impose independent constraints for understanding the nature 
of this "dark particle accelerator" . 
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1. INTRODUCTION 

Recent surveys of the central region of our 
Galaxy with the High Energy Steroscopic System 
(H.E.S.S.) have uncovered a number of 7— ray 
sources in the TeV regime (Aharonian et al. 2002, 
2005a, 2005b, 2006a). Different from the cases 
of pulsar wind nebulae (PWNe) and young super- 
nova remnants (SNRs), some of these sources have 
no non-thermal X-ray counterpart yet been iden- 
tified. Among them, HESS J1745-303 is one of the 
most enigmatic objects. 

HESS J1745-303 was firstly discovered by the 
H.E.S.S. Galactic Plane Survey (Aharonian et al. 
2006a) and was subsequently investigated in de- 
tails with dedicated follow-up observations (Aha- 
ronian et al. 2008). The TeV 7— ray image shows 
that it consists of three spatial components (i.e. 
Regions A, B and C in Fig. 1 of Aharonian et 
al. 2008). Owing to the lack of spectral variabil- 
ity and the insignificant dip among these regions 
in the existing data, it was argued that they are 
originated from a single object (Aharonian et al. 
2008). This inference suggests that HESS J1745- 



303 as one of the largest unidentified TeV sources 
which has an angular size of ^ 0.3° x 0.5°. 

Searches for the possible non-thermal diffuse 
X-ray component in region A of HESS J 1745- 
303 have been conveyed with XMM-Newton and 
Suzaku (Aharonian et al. 2008; Bamba et al. 
2009) . None of these observations have resulted in 
any evidence for the diffuse X-ray emission. This 
imposes a TeV-to-X-ray flux ratio larger than ^ 4 
(Bamba et al. 2009) which is larger than the typi- 
cal value of PWNe and SNRs (i.e. less than 2) (cf. 
Matsumoto et al. 2007; Bamba et al. 2007, 2009). 
Because of the non-detection of counterparts in X- 
ray/radio, HESS J1745-303 is dubbed as a "dark 
accelerator" (Bamba et al. 2009). 

While no non-thermal diffuse X-ray emission 
has yet been found, a possible excess of neutral 
iron line emission was discovered in the direction 
toward the region A of HESS J 1745-303 (Bamba 
et al. 2009). Together with its proximity to the 
Galactic center and the positional conicidence of 
a molecular cloud (Aharonian et al. 2008), the 
line emission is suggested to be the reflected X- 
rays originated from the previous activity in the 
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Galactic center (Bamba et al. 2009). This molec- 
ular cloud can be interacted with the shock from a 
nearby SNR G359.1-0.5 (Bamba et al. 2000, 2009; 
Lazendic et al. 2002; Ohnishi et al. 2011) and 
produce the observed 7— rays through the acceler- 
ation of protons and/or leptons (see Bamba et al. 
2009). Nevertheless, this proposed scenario can- 
not be confirmed unambiguously. In view of the 
presence of many surrounding objects (see Fig. 1 
in Aharonian et al. 2008), including the "mouse" 
pulsar (i.e. PSR J174 7-2958), one cannot rule out 
these objects as the source of energetic particles 
simply based on the TeV results (Aharonian et al. 
2008). Furthermore, with the current information, 
it is not possible to discriminate hadronic model 
and Icptonic model (see Fig. 5 & 6 in Bamba et 
al. 2009). In order to do so, investigations in lower 
energy regime are required. 

It is interesting to note that HESS J 1745- 
303 is positionally coincident with an unidenti- 
fied EGRET source 3EG J1744-3011 (Hartman 
et al. 1999). Different from HESS J1745-303, 
3EG J1744-3011 was suggested to demonstrate 
long-term variability (Torres et al. 2001). Also, 
based on the MeV-GeV spectrum observed by 
EGRET, the extrapolated fiux in the TeV regime 
overshoots that observed by H.E.S.S.. There- 
fore, Aharonian et al. (2008) considered that 
3EG J1744-3011 is unrelated to HESS J1745-303. 

After the commence of the Large Area Tele- 
scope (LAT) onboard Fermi Gamma-ray Space 
Telescope, a detailed investigation of this dark ac- 
celerator in the MeV— GeV regime is now feasible 
with its much improved spatial resolution and sen- 
sitivity. However, among 1451 objects detected 
by LAT during the first 11 months, we do not 
identify any source corresponding to HESS J 1745- 
303 / 3EG J1744-3011 (Abdo et al. 2010). Very 
recently, in an analysis of the 7— rays from the 
Galactic center with first 25 months LAT data, 
a new serendipitous source was found to coincide 
spatially with HESS J1745-303 / 3EG J1744-3011 
(Chernyakova et al. 2011). In this paper, we re- 
port a detailed analysis of this source with LAT 
observation in the first ~ 29 months. 

2. DATA ANALYSIS & RESULTS 

In this analysis, we used the data obtained by 
LAT between 2008 August 4 and 2010 Decem- 



ber 23. The Fermi Science Tools v9rl8p6 pack- 
age is used to reduce and analyze the data in the 
vicinity of HESS J1745-303. Only the events that 
are classified as class 3 or class 4 are adopted. 
The post-launch instrument response functions 
"P6_V3_DIFFUSE" were used throughout this in- 
vestigation. 

With the aid of the task gtlike, we performed 
unbinned maximum-likelihood analysis for a cir- 
cular region-of-interest (ROI) with a 10° diameter 
centered on the nominal position of HESS J 1745- 
303 (i.e. RA=17''45'"1.999^ Dec=-30°22'l2.0" 
(J2000)). The size of ROI have been chosen to 
avoid the surrounding bright sources so as to re- 
duce the systematic uncertainties due to the in- 
accurate background subtraction in this complex 
region. For subtracting the background contri- 
bution, we included the Galactic diffuse model 
(gll_iem_v02.fit) and the isotropic background 
(isotropic_iem_v02.txt), as well as 41 sources in 
the first Fermi /LAT catalog (IFGL; Abdo et al. 
2010) within 10° from the aforementioned center. 

To begin with, we compared the spectral prop- 
erties inferred by LAT and EGRET. For a con- 
sistent comparison with Hartman et al. (1999), 
we used events with energies > 100 MeV for our 
initial analysis. We assumed a power-law (PL) 
spectrum for HESS J1745-303 as weh as all IFGL 
sources in our consideration. All the sources are 
assumed to be point sources throughout this inves- 
tigation. The best-fit model yields a photon index 
of r = 2.16 ± 0.0^ and a test-statistic (TS) value 
of 499 which corresponds to a significance of 22(7. 
This is consistent with the significance reported 
in the preliminary analysis by Chernyakova et al. 
(2011). The photon index is found to be consistent 
with 3EG J1744-3011 (i.e. P = 2.17 ±0.08). To 
further compare with the EGRET results, we com- 
pute the integrated photon flux in 0.1 — 10 GeV 
which is found to be 2.06t!^;23 x 10"^ cm^^ s~i. 
This is 3 times smaller than that of 3EG J 1744- 
3011 (i.e. 6.39 ± 0.71 x lO"'^ cm'^ s'^). The 
discrepancy can be due to the improved spatial 
resolution of LAT, and hence the estimation of 
background contribution from the nearby sources 
is more accurate than EGRET. 

Considering the background, as HESS J1745- 



^AU errors quoted in this paper are statistical only and are 
computed for a confidence interval of la. 
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303 is located close to the Galactic center where 
the diffuse 7— ray emission is very intense, the con- 
tamination can possibly be large at lower energies. 
Also, the point spread function (PSF) of LAT is 
narrower at higher energies. While the 68% con- 
tainment radius at 100 MeV is ~ 4.5°, it is only 
~ 0.8° at 1 GeV. E Therefore, the contamina- 
tion due to the PSF wings of the nearby sources 
can also be minimized by limiting the analysis at 
higher energies. In order to minimize the system- 
atic uncertainty in the background modeling so 
as to obtain robust results, we restricted the sub- 
sequent analysis in 1 — 20 GeV. In this adopted 
band, the best-fit PL model results in a TS value of 
332 which corresponds to a detection significance 
to be ~ 18(T. To examine the robustness of the 
detection, we have also considered the systematic 
uncertainty of the Galactic diffuse emission back- 
ground. Following Abdo et al. (2010), we repeat 
the analysis by varying the slope and the normal- 
ization of the Galactic diffuse model in — 0.07 
and ±10% respectively. Within the uncertainty of 
the background model, the detection significance 
remains over lOcr. This best-fit PL model yields a 
photon index of F = 2.60 ±0.05 and an energy fiux 

1 in 1-20 GeV. 
-ray spectrum as 
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Fig. 1.— Fermi LAT spectrum of HESS J1745- 
303. The solid line represents the best- fit power- 
law model. 

perturbations in the parameter space. In view of 
the problem of the convergence, we will not further 
consider BKPL model in all subsequent analysis. 

Although the extra parameter, -Bcutoff, in the 
PLE model is not statistically required, one is not 
able to completely rule out it in view of its capa- 
bility in depicting the observed data. As a PLE 
model typically describes the 7— ray spectrum of 
a pulsar (Abdo et al. 2010b, 2011), we speculate 
if there is any hidden pulsar in HESS J1745-303. 
To test this hypothesis, we have performed a blind 
search for coherent pulsation. The arrival times of 
each photon were barycentric corrected with the 
nominal position of HESS J1745-303 adopted. To 
minimize the impact due to the ignorance of the 
spin-down rate, we have divided the full time span 
of the adopted data into 5 segments and run the 
Fourier analysis in each segment independently 
with the aid of the tool gtpspec. From each com- 
puted power spectrum, we picked out 10 peaks 
and investigated if there is any correlation among 
different segments. However, we did not identify 
any promising periodicity candidate in this data. 
Hence, we conclude that there is no evidence for a 
hidden pulsar in HESS J1745-303. 

As 3EG J1744-3011 was reported to be a vari- 
able (Torres et al. 2001), we also examine the 
variability with LAT data. First, we extracted 
the light curve obtained from the data within 1° 

glast.slac.stanford.edu/software/IS/glast_lat-performanfeCtaiithe nominal position of HESS J 1745-303 with 



of 5.25j:}:^^ X 10"" ergs cm-^ s" 
The best-fit PL model and the 7 
seen by Fermi LAT is shown in Figure [1] 

Besides a simple PL model, we have also exam- 
ined if a broken power-law (BKPL) or an expo- 
nential cutoff power-law (PLE) can describe the 
spectrum better. The fittings with BKPL and 
PLE results in the TS values of 335 and 333 re- 
spectively. Based on the likelihood ratio test, the 
additional spectral parameters in BKPL/PLE are 
not strongly required which suggests that a sin- 
gle PL is statistically sufficient to describe the 
data. For the PLE model, the best-fit photon in- 
dex and cut-off energy are F = 2.31 ± 0.07 and 
i;^^t„ff = 12.63 ± 4.69 GeV respectively. On the 
other hand, for the BKPL, the initial fitting re- 
sulted in a break energy and the photon indicies 
of iJbrcak = 3.68 ± 0.22 GeV, Fi = 2.32 ± 0.15 and 
F2 = 3.09 ± 0.26 respectively. However, different 
from the best-fit solutions inferred from PL and 
PLE models, we found that the solution inferred 
from the BKPL fit is unstable subjecting to the 
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Fig. 2. — (a) Test-statistic (TS) map in 1 — 20 GeV of a region of 2° x 2° centered at the nominal position 
of HESS J1745-303. The color scale that used to indicate the TS value is shown by the scale bar below. 
The blue circle represents the la positional error circle determined by gtfindsrc. Various TeV emission 
components of HESS J1745-303 (i.e. regions A, B and C in Fig. 1 of Aharonian et al. 2008) are illustrated 
by the black dashed circles, (b) Same as Fig. 2a but in the energy range of 10 — 20 GeV. 



a binning factor of 10 days in the energy range of 
1 — 20 GeV. By fitting a horizontal line to the light 
curve, we obtain a reduced chi-square of = 1.5 
for 87 degrees of freedom. Hence, there is no 
strong evidence for any variability or flaring. 

For a further investigation of the possible spec- 
tral and flux variability, we divided the whole data 
into flve segments of equal time span and per- 
formed an unbinned likelihood analysis on each 
segment. The results are summarized in Table [T] 
A simple PL was adopted for all the fittings. 
Within the tolerance of the statistical uncertain- 
ties, we conclude that neither the spectral shape 
nor the flux varies among these segments. 

We have computed the 2° x 2° TS map in 
1 — 20 GeV centered at the nominal position of 
HESS J1745-303 by using the tool gttsmap. This 
is displayed in Figure [2^. Utilizing gtfindsrc, we 
determined the best-fit postion in 1 — 20 GeV to be 
RA=17''44™31.440" Dec=-30°20'32.28" (J2000) 
with a 1(7 error radius of 0.05°. The position lo- 
cates between regions A and C of HESS J1745-303. 
Given that the PSF has a 68% containment radius 



of ^ 0.8° at 1 GeV, the source extent inferred in 
this band is consistent with that of a point source. 
The relatively wide PSF do not allow us to de- 
termine whether GeV emission is associated with 
any particular TeV feature. To further examine 
the spatial nature, we also computed the TS map 
in 10 — 20 GeV which is displayed in Figure 
Since the 68% containment radius at 10 GeV is 
~ 4 times smaller than that at 1 GeV, the feature 
can be better resolved. We found that the peak TS 
value found in this band is 41 which corresponds 
to a significance of ~ 6(j. Within the systematic 
uncertainty of the Galactic diffuse background, we 
found that the detection significance of the source 
remains over Aa in this band. The best-fit position 
in 10 - 20 GeV is found to be RA=17''43™44.160" 
Dec=-30°26'24.00" (J2000) with a la error ra- 
dius of 0.05°. This differs from that inferred in 
1 - 20 GeV by 0.2°. We note that the GeV fea- 
ture found in this hard band is apparently peaked 
at region C and possibly extended to region B. Al- 
though it appears to be extended in the hard band, 
the relatively low detection significance in this en- 
ergy range does not allow a firm conclusion. 
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Table 1 

7-ray spectral properties and energy fluxes of hess j1745-303 at different epochs. 



Time segment 


r 


f^%l - 20 GcV) 


Mission elapsed time (s) 




10^^^ ergs cm^2 g-i 



239557417- 


-254601478.2 


2.67 ±0.23 


4.45i^^,^^ 


254601478.2- 


-269645539.4 


2.45 ±0.18 


6.65;^:'3f 


269645539.4- 


-284689600.6 


2.57 ±0.05 


5.73+ff 


284689600.6- 


-299733661.8 


2.70 ±0.25 




299733661.5 


i-314777723 


2.80 ±0.30 


4 oc+ig.ie 



^The quoted errors of energy have taken the statistical uncer- 
tainties of both photon index and prefactor into account. 



3. DISCUSSION 

In this paper, we have reported a detailed study 
of HESS J 1745-303 with Fermi LAT data which 
provides a missing piece in understanding the na- 
ture of this dark accelerator. In view of the pu- 
tative variability of 3EG J1744-3011, Aharonian 
et al. (2008) argued that the GeV counterpart 
is unlikely to be associated with HESS J 1745- 
303. However, we do not find any evidence for 
the spectral/flux variability from the LAT data. 
The discrepancy of the LAT result and that in- 
ferred from EGRET is most likely due to the dif- 
ferences between their instrumental performance. 
With the improved angular resolution and sensi- 
tivity of LAT, many previously unknown sources 
in the proximity of HESS J1745-303 are now de- 
tected. As several sources within 2° are found to 
be variable, including IFGL J1747.2-2958 which is 
the 7-ray counterpart of PSR J1747-2958 (Abdo 
et al. 2010), the source confusion in the EGRET 
data can possibly lead to the apparent variability. 

It is interesting to compare the spectral prop- 
erties inferred in the GeV regime with those ob- 
tained by H.E.S.S.. First, the photon index in- 
ferred from the LAT data (i.e. F = 2.60 ± 0.05) 
is similar to that inferred in TeV (cf. Tab. 2 of 
Aharonian et al. 2008). Furthermore, the extrapo- 
lated photon flux in 1 - 10 TeV with the best-fit PL 
model is found to be (2.21 1^;^;^) x IQ-^^ cm-^ s'^ 
which is consistent with any individual spatial 
component observed by H.E.S.S. within la uncer- 
tainties (cf. Tab. 2 of Aharonian et al. 2008). 

Since the TeV spectral data of region A is avail- 



able (cf. Fig. 2 of Aharonian et al. 2008) Q, we fur- 
ther compare this particular region with the GeV 
spectrum by constructing a spectral energy distri- 
bution (SED) which is display in Figured Both 
data can be simultaneously fitted with a single PL 
of r = 2.63 ±0.03. This clearly demonstrates that 
the TeV spectrum of this spatial component can 
be smoothly connected to the GeV spectrum. 

For regions B and C, despite that no TeV 
spectral data is currently available, we note that 
their spectral shapes reported by Aharonian et 
al. (2008) are similar to that of region A. Within 
the statistical uncertainties of the spectral prop- 
erties inferred in both GeV and TeV regimes, the 
TeV spectra of these regions can also be possi- 
bly connected to that of GeV. Further investiga- 
tion by H.E.S.S. is strongly encouraged, particu- 
larly for region C as most of the 7— ray emission in 
10 — 20 GeV apparently coincides with this spatial 
component. 

As both PL and PLE model can describe the 
GeV spectrum equally well, we cannot discrimi- 
nate these competing models based on the current 
data. In view of the exponential cut-off, the spec- 
tral connection between the GeV and TeV regimes 
cannot be established with the PLE model. How- 
ever, we would like to point out that the best-fit 
cut-off energy, E^utoS = 12.63 ±4.69 GeV, falls in 
the highest energy bin of the LAT spectrum (cf. 
Fig. [T] and Fig. [3|). Owing to the small photon 
statistics in the hard band, the statistical uncer- 
tainty of the highest energy bin is rather large. 



^http: / /www. mpi-hd.mpg.de/hfm/HESS 
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Fig. 3. — Spectral energy distribution of 
HESS J1745-303 as observed by Fermi LAT and 
H.E.S.S. The H.E.S.S. spectrum is for region A. 
The soUd hne represents the best-fit power-law 
model inferred in the joint analysis of both data 
sets. 

In view of this, it remains unclear whether the 
inferred cut-off is genuine as this particular data 
point is sensitive to the systematic uncertainty of 
the background. Analysis of the LAT data in 
higher energies (e.g. > 10 GeV) with sufficient 
photon statistic in the future can help to discrim- 
inate these two models. 

The detection of 7— rays provides a strong ev- 
idence for the particle acceleration. It should 
be noted that there are two SNRs and two 
pulsars around HESS J1745-303 which can be 
the potential high energy particle injector (cf. 
Fig. 1 of Aharonian et al. 2008). Based on the 
small distances of G359.0-0.9 and PSR B1742- 
30, Bamba et al. (2009) argued that these two 
sources are foreground objects which are unre- 
lated to HESS J1745-303. With the GeV coun- 
terpart revealed by LAT, we can now safely 
rule out the possibility that PSR B1742-30 is 
the major contributor as its spin-down flux, 
E/ATidP ^ 10~^^ erg cm~^ s~^, is lower than 
the energy flux observed in the GeV regime (cf. 
Manchester et al. 2005). On the other hand, 
for PSR J1747-2958, the sum of 7-ray flux of 
IFGL J1747.2-2958 and HESS J1745-303 observed 
by LAT only consumes ^ 5% of its spin-down 
flux ('^ 3 X 10~^ erg cm~^ s^^). Therefore, it 
is energetically possible to be the source for the 
high energy particles. Nevertheless, with its east- 



ward proper-motion (Hales et al. 2009), its back- 
ward extrapolated position by its spin-down age 
is found to have an offset of ~ 0.6° and ^ 0.8° 
from the best-fit positions inferred in 1 — 20 GeV 
and 10 — 20 GeV respectively. If the feature found 
by H.E.S.S. and LAT is indeed a PWN, then it 
is one of the most peculiar system because of its 
large positional offset with respect to the pulsar 
position. 

While the contribution of PSR J1747-2958 is 
uncertain, the interaction between the shock from 
G359. 1-0.5 and a molecular cloud in its neighbour- 
hood is considered to be a more viable means to 
produce the observed 7— rays. For both hadronic 
and leptonic scenarios, Bamba et al. (2009) have 
modeled the broadband spectrum of the region A 
of HESS J1745-303. However, comparing Fig. [1] in 
this paper and Fig. 5 & 6 in Bamba et al. (2009), 
the GeV flux predicted in all the scenarios consid- 
ered in their work are at least an order of mag- 
nitude lower than that observed by LAT. Revised 
theoretical investigation with constraints provided 
by LAT is therefore required. Moreover, while 
the other known systems of SNRs interacting with 
molecular clouds have their remnant shells found 
in other wavelengths coincided with the 7— ray 
emission (e.g. Castro & Slane 2010; Abdo et al. 
2010c), there is no such evidence for HESS J1745- 
303. We should point out that all the previous 
multiwavelength campaign were targeted only at 
the bright component of HESS J1745-303, namely 
the region A. Since the LAT observation suggested 
that the GcV emission can possibly be originated 
from the regions C/B, observational investigations 
aim at these regions are encouraged for a further 
understanding of this mysterious object. 

Based on the above discussion, we have to ad- 
mit that the energy injection source of HESS J1745- 
303 remains unclear. It is instructive to compare 
it with other nearby high energy sources. Aha- 
ronian et al. (2006b) have reported observations 
of an extended region of very high energy (VHE, 
> 100 GeV) 7-ray emission correlated spatially 
with a complex of giant molecular clouds in the 
central 200 pc of the Milky Way. It appears that 
TeV emissions from the molecular clouds in the 
vicinity of the Galactic Center are quite common 
phenomena. In addition, similar to the case of 
HESS J1745-303, 6.4 keV lines have been com- 
monly detected from many molecular clouds near 
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the Galactic Center, e.g. Sgr B2 (Koyania et al. 
1996; Murakami et al. 2000), Sgr C (Murakami 
et al. 2001) and others (Bamba et al. 2002; Pre- 
dehl et al. 2003; Nobukawa et al. 2008; Nakajima 
et al. 2009; Koyama et al. 2009). It has been 
speculated that these neutral iron lines arise from 
the reflection by the dense molecular clouds which 
are irradiated by the nearby X-ray sources (e.g. 
Koyama et al. 2007; Inui et al. 2009). On the 
other hand, Dogiel et al. (2009) argued that these 
6.4 keV lines from the molecular clouds are ex- 
cited by a background intensity of subrelativistic 
protons coming from the escaped part of a past 
captured star by the Galactic supermassive black 
hole Sgr A*. The periodic stellar capture events 
may explain the recent observed Fermi bubble 
(Cheng et al. 2011a). This subrelativistic proton 
wind can also form shock by hitting the clouds and 
produce relativistic protons. It is shown that the 
decay of neutral pions produced by hadronic colli- 
sions between the accelerated relativistic protons 
in the clouds can emit a power-law 7— ray spec- 
trum from 30 MeV to 10 TeV without any spectral 
break (Cheng et al. 2011b in preparation). If this 
is true, then the high energy emission from vari- 
ous molecular clouds in the vicinity of the Galactic 
center are correlated to the past activities of the 
Galactic center and their intensity might be corre- 
lated with the propagation history of the injection 
of these subrelativistic protons escaped from the 
capture past events. 
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